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ABSTRACT 

The progenitor of SN 2005cs, in the galaxy M51, is identified in pre-explosion HST ACS 
WFC imaging. Differential astrometry, with post-explosion ACS HRC F555W images, per- 
mitted the identification of the progenitor with an accuracy of 0.006". The progenitor was 
detected in the F814W pre-explosion image with / = 23.3 ± 0.2, but was below the detection 
thresholds of the F435W and F555W images, with B < 24.8 and F < 25 at Sct. Limits were 
also placed on the U and R band fluxes of the progenitor from pre-explosion HST WFPC2 
F336W and F675W images. Deep images in the infra-red from NIRI on the Gemini-North 
telescope were taken 2 months prior to explosion, but the progenitor is not clearly detected on 
these. The upper limits for the JHK magnitudes of the progenitor were J < 21.9, H < 21.1 
and K < 20.7. Despite having a detection in only one band, a restrictive spectral energy dis- 
tribution of the progenitor star can be constructed and a robust case is made that the progenitor 
was a red supergiant with spectral type between mid-K to late-M. The spectral energy distri- 
bution allows a region in the theoretical HR diagram to be determined which must contain the 
progenitor star The initial mass of the star is constrained to be Mzams = ^-2^&' which is 
very similar to the identified progenitor of the type II-P SN 2003gd, and also consistent with 
upper mass limits placed on five other similar SNe. The upper limit in the deep K-band image 
is significant in that it allows us to rule out the possibility that the progenitor was a signifi- 
cantly higher mass object enshrouded in a dust cocoon before core-collapse. This is further 
evidence that the trend for type II-P SNe to arise in low to moderate mass red supergiants is 
real. 
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1 INTRODUCTION 

Stellar evolution models predict that stars with initial masses > 
8Mq end their lives in a core-collapse induced supernova (CCSN), 
and the CCSN mechanism is also related to the energetic gamma- 
ray burst (GRB) phenomenon ( Matheson et al. 2003). The study 
of the progenitors of CCSNe is essential to both the understand- 
ing of the evolution of SNe and the massive stars which produce 
them. Fortuitous pre-explosion observations of the progenitor can 
allow the determination of the evolutionary state of the star prior 
to SN explosion. The progenit or of SN 1987A was found to be 
a blue supergiant ( Walbom et all 198^ . whereas the progenitor of 
SN 1993 J was identified as a red sup ergiant with a UV excess due t o 
a hot binary companion ( Alderina e t all 1994 iMaund et alliooj) . 
In the cases of SNe 1987 A and 19931 the direct observations of the 
progenitors have helped explain peculiarities in the SN explosions 
themselves. The progenitor of the type II-P SN 2003gd was iden- 
tified as a red supergiant, the canonical prediction of stellar ev olu- 
tion m o dels for t his typ e of SN ( Smartt et al. 2004; Van Dvk e tal] 
l2003lJ) . lLi et'ZI<2005ii) report the identification of a yellow super- 



giant proge nitor for the type II - P SN 2004et. Recentlv |Maund et alj 
( 2005) and iMaund & SmarttI <2005h have placed limits on CC- 
SNe, of various types, whose progenitors were not detected on pre- 
explosion imaging. Maund & Smartt (2005) demonstrate the im- 
portance of high resolution post-explosion imaging of SNe, using 
differential astrometry to locate (with sub-pixel accuracy) the pro- 
genitor on pre-explosion images. 

In this letter we report the identification of the progenitor 
of the Type II-P SN 2005cs in M51, the "Whirlpool Galaxy." 
SN 2005cs was discovered by Kloeht (2005) on 2005 June 
28.9. The SN was spectroscopically classified as a Type II by 
iModiaz et al! J2w3), with P-Cygni profiles o f the B aimer and He 
lines.lkichmond ( 20o|) and lLi et all J2005d) reported the identifi- 
cation of two different stars as possible progenitors for this object: a 
blue and a red star respectively. Li et al. ( 2005b) give the position 
of SN 2005CS as aaooo = 13''29"'52.=76,52000 = 47°10'36."ll. 
This implies an offset between SN 2005cs and the centre of M5 1 
of Aa = 6" and A J = -64.7", or 2.7kpc at a distance of 8.4Mpc 
Jpeldmeier et alll99l with PA. and inclination of M5 1 as quoted 
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Table 1. Pre-and post-explosion observations of the site of SN 2005cs used 
in this paper. 



Date 


Instrument 


Filter 


T-* Y n T'l m p 

(s) 




Pre-explosion images 




2005 Jan 20-21 


ACSAVFC 


F435W 


2720 




ACSAVFC 


F555W 


1360 




ACSAVFC 


F814W 


1360 


1999 Jul 21 


WFPC2 


F336W 


1200 




WFPC2 


F555W 


1200 




WFPC2 


F675W 


500 


2005 Apr 29 


NlRl 


J (G0202) 


500 




NlRl 


H (G0203) 


500 




NIRl 


K (G0204) 


625 


2005 Jul 24 


Post-explosion images 

ACS/HRC F555W 


80/160/1704 



by HyperLEDA^ as 163° and 47.5° respectively). While this paper 
was in the process of submission, a similar study was submitted for 
publication by Li et al. ( 2005c). There are some differences in the 
results of the two studies and we briefly highlight them in this letter 
where appropriate. 



2 OBSERVATIONS AND DATA REDUCTION 

The pre- and post-explosion images analysed here are summarised 
in Table Q Pre-explosion images of the site of SN 2005cs were 
available in the Hubble Space Telescope (HST) archive^, with 
multi-colour imaging with the Wide Field Planetary Camera 2 
(WFPC2) and Advanced Camera for Surveys (ACS) Wide Field 
Channel (WFC). Individual HST WFPC2 observations were re- 
trieved from the STScI archive. These observations were calibrated 
with the On-the-fly-recalibration (OTFR) pipeline. The site of SN 
2005CS was located on the WF2 chip in the 1999 WFPC2 obser- 
vations, which has 0.1" pix~^. The WFPC2 observations were 
combined for the extraction of cosmic rays and corrected for ge- 
ometric distortions. Aperture photometry was then conducted on 
these images using DAOphot. PSF photometry was cond ucted on 
the W FPC2 observations using the HSTPHOT package (Dolphin 
I2OOQ) . HSTphot includes corrections for chip-to-chip variations, 
charge transfer efficiency and transformations from the WFPC2 
photometric system to the standard Johnson-Cousins magnitude 
system. The WFPC2 pre-explosion observation complemented the 
ACS pre-explosion imaging with observations in F336W and 
F675W (approximately U and R respectively). 

The ACS pre-explosion images were acquired as part of pro- 
gram GO-10452 (PI: S.V.W. Beckwith), which conducted a large 
deep mosaic of M51 (NGC 5194) and its companion NGC 5195. 
These observations were in four bands F435W, F555W, F814W and 
F658N . The drizzled combin ed mosaic frames were available from 
STScI jMutchler et alj2005h .with scale 0.05" pix"\ Aperture and 
PSF photometry was conducted on the mosaic frames using the 
PyRAF implementation of the DAOPHOT algorithm. Empirical 



^ http : / /www — obs.univ — lyonl.fr/hypercat/ 
^ http://archive.stsci.edu/hst/ 



aperture corrections were applied to the data, with the charge trans- 
fer efficiency corrections of iRiess & Mack ( 20oJ) and the colour 
transformation equations and updated zeropoints of ISirianni et alJ 
(2005). 

Infrared pre-explosion observations of the site of SN 2005cs 
were conducted with the Gemini North Telescope Near Infra- 
red Imager (NIRI) in the JHK bands. These observations were 
conducted as part of program GN-2005A-Q-49 (PI: S. Smartt). 
The f/6 camera of NIRI was used to provide imaging with scale 
0.117" pix^^. The observations were split into a series of short 
sub-exposures, to overcome the sky background, with small offsets 
between each sub-exposure. On completion of the series of obser- 
vations on target, a comparable series of short exposures at a large 
offset was immediately acquired to measure the levels of the IR sky 
background. The NIRI images were reduced and combined, with 
sky subtraction, using the IRAF gemini package of reduction tools 
for the NIRI instrument. Photometry on the JHK images was done 
by PSF fitting within the IRAF implementation of DAOPHOT. The 
zero-point of the NIR filters was determined using two standard 
stars (UKIRT Faint Standards 133 and 136) taken before and after 
the science frames. They gave comparable results, and suggested 
that the period of the observations was stable and photometric. The 
image quality was 0.5" in K and 0.6" in J and H. 

Post-explosion ACS High Resolution Channel (HRC) imaging 
was conducted as part of program GO- 10498 (PI: S. Smartt). These 
observations were only acquired in the F555W, as good colour in- 
formation on the surrounding stellar population was possible from 
the pre-explosion observations. These were photometered in a sim- 
ilar way to the pre-explosion ACS/WFC observation using PyRAF 
DAOPHOT. The same corrections were applied, although the mag- 
nitudes were left in the ACS VegaMag photometric system. 

The geometric transformation between the F555W images of 
WFPC2 1999 observations with 2005 ACS/WFC observations was 
calculated, using the iRAF task geomap, with 1 1 common stars. X- 
and Y- offsets between the F555W image and other filter images 
from the same observing run were calculated by cross-correlating 
the images, using the Stsdas crossdriz task with the offset calcu- 
lated from the cross correlation image using shiftfind. The geomet- 
ric transformation was also calculated between the 2005 ACS/WFC 
F814W image and the NIRI J-band image, with the shift calcu- 
lated between the ACS/WFC F555W and F814W images applied 
to bring the transformed NIRI J-band image to the same coordi- 
nate system as the ACS/WFC F555W image. Similarly to the HST 
pre-explosion observations, offsets between the J-band image and 
the H and K images were calculated using the cross correlation 
method. The uncertainty in transforming the position of a star from 
one image to another was calculated as the sum in quadrature of 
the positional uncertainties in the reference image and transformed 
image and the r.m.s. uncertainty of the transformation, by geomap. 
The position uncertainties on the individual images were assessed 
by measuring the scatter in the positions of the stars, measured us- 
ing the four centering algorithms of DAOPHOT, of all the stars 
used to calculate the geometric transformation. The individual im- 
age position uncertainty was taken as the sum in quadrature of the 
mean scatter and standard deviation about this mean in both the X 
and Y directions. 

3 OBSERVATIONAL RESULTS 

SN 2005cs was detected in the post-explosion ACS/HRC F555W 
imaging with mp555w = 15.28 ± 0.01. The position of SN 2005cs 
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on the post-explosion frame was calculated to within ±0.002" . The 
position of the SN on pre- and post-explosion images in shown on 
Fig. 1. The SN position on the pre-explosion ACSAVFC F555W 
image was identified to within 0.006" (which is the combined 
uncertainties on the determination of positions on the pre- and 
post-explosion images and the geometric transformation). A pro- 
genitor star is significantly detected on the F814W image with 
niF8i4w = 23.26 ± 0.03. There is no object in the pre-explosion 
F435W and F555W images, and we estimate 5(t detection lim- 
its of mF435w = 24.9 and mpsssw = 25 (using the methods of 
Maund & Smartt 2005). In order to add further limits to the spec- 
tral energy distribution (SED) of the progenitor, we determined 
the limiting magnitude of the WFPC2 F336W and F675W im- 
ages. Due to the lower resolution of these images compared to 
the ACS, the progenitor position was too close to the brighter ob- 
ject to the southwest to allow accurate simultaneous PSF fitting. 
Hence a large (6 pixel) aperture flux was determined at the pro- 
genitor position, and the flux of stars which had well determined 
individual small aperture photometry was subtracted. This sug- 
gested there was no residual flux fro m^ any object brig hter than 
niF336W = 21.1 and mpersw =23.6. Li et al. (2005c) report a 
magnitude of J = 24.15 ± 0.20 for the progenitor which is signifi- 
cantly fainter than our measurement. We have extensively checked 
our ACS photometry and zeropoints using several methods (e.g. 
DOLphot"^ and several alternative experiments with PSF and aper- 
ture photometry) and also compared the F814W magnitudes (from 
HSTphot and our own estimates) of faint isolated stars in WFPC2 
archive images with our ACS measured magnitudes. We find no 
problems or inconsistencies in our analysis and confirm the magni- 
tude as mF8i4w = 23.26 ± 0.03. 

The location of the SN on the pre-explosion Gemini NIRI 
imaging was identified to within 0.078". There is a JHK source 
clearly detected close to the position of the J-band progenitor, how- 
ever its position is more than 2g away and hence is not spatially 
co incident with the progenitor. This object has also been detected 
bv lLietaljEoOSd) in NICMOS pre-explosion images, and is coin- 
cident with the 7-band star to the north west of the progenitor. At 
the position of the progenitor in the JHK NIRI images, the flux 
does appear somewhat higher than the surrounding background, but 
given the ground-based resolution there is no evidence for detection 
of a single source (see Fig. I). A further check was performed to de- 
termine that no object could be detected in the 7^-band. The bright 
objects on the image were subtracted off by PSF fitting, and fainter 
objects were identified on the resultant frame. After three iterations 
of PSF fitting and subtraction, a boxcar median filtered image was 
constructed, which should be representative of the varying galaxy 
background. This was subtracted off the original frame and the pho- 
tometric procedure was repeated. There was no evidence of a single 
point source at the progenitor position. The detection limits in the 
J H K-b&nA images where determined at the SN position. In the 
first instance statistical calculations of the 3(7 detection limits were 
determined from the detector characteristics and the sky noise. This 
was compared with the measured magnitudes of the faintest real 
objects in the images which could be visually identified and which 
were both well fit with a PSF and subtracted off cleanly. The 3(7 cal- 
culated limits were significantly fainter (by ~0.7 mag) than these 
objects, but the 5(7 limits agreed well with the magnitudes of the 
faintest sources measured in the image. Hence we adopted the 5(7 
limits, which are: J = 21.9, H = 21.1 and K = 20.7. The JH 



^ http://purcell.as.arizona.edu/dolphot/ 



limits are consistent with the NICMOS limits reported bv lLi et all 
( 2005c), although the additional A"-band upper limit we present has 
significant implications (see Section 4.3). 



4 DISCUSSION 

4.1 Estimates of metallicity of the progenitor region 

HII regions around M5I were spectroscopically studied by 
iBresolin et alj i2004h for the purposes of abundance analysis. The 
sight line corrected angular distance of SN 2005 cs is 66.5" or 
0.21i?o (where Rq is the characteristic radius which lsresolin et all 
(2004) give as 5.4'). The oxygen abundance, at the radius of 
SN 2005CS, was calculated to be 12 + log(0/H) = 8.66 ± 0.11, 
which is similar to abundance in the Hll region CCM56 
(>8.58 dex) closest to SN 2005cs. Eresohn et al. ( 2004) point out 
that their new abundance determinations, which are based on di- 
rectly determined electron temperatures, are significantly lower 
than previous estimates in M51 and that the previously reported 
very metal rich abundances are likely to have been inaccurate. As 
^SDlund et al. (2004) give the solar oxygen abundance as 8.66 ± 
0.05, this implies that the metallicity appropriate for SN 2005cs 
and its progenitor star is very likely to be in the range 1 ± 0.3^0 



4.2 Estimates of Reddening 

iModiaz et alj i2005h reported the presence of Galactic and host 
galaxy components of the NalD lines in an early spectrum of SN 
2005cs. They reported that the equivalent width of the host-galaxy 
compone nt was similar in s trength to the Galactic component 
of O.2A. iTuratto et alJ 1200 3^ provide relationships between the 
strength of the NalD absorption and the reddening from the inter- 
vening medium. If one assumes that the strengths of the two NalD 
components may be added then the total NalD line strength implies 
a red dening between E{B — V) — 0.05 — 0.16. Bresolin et ^ 
(2004) record the reddening coefficient c(H/3) for the nearby HII 
region CCM 56 of 0.23 which c orresponds to E{B - V) = 0.16 
with a Galactic iHowarthI il983i) reddening law. The three colour 
photometry of red supergiants (B — V > 0.2) within 2" of SN 
2005cs was also used to estimate the reddening as detailed in 
iMaund & Smartll i2005t) . Stellar positions were plotted on a two- 
colour diagram {B — V ys.V — I) and compared with a theoretical 
supergiant colour sequence. The weighted mean displacement of 
the stars' positions, along the reddening vector, from the theoretical 
colour sequence provided a measure of reddening. This technique 
yielded a value of E{B — V) = 0.12 ± 0.01, and we hence adopt a 
reddening of E{B - V) = 0.14 ± 0.02 throughout the rest of this 
paper. 



4.3 The Progenitor of SN 2005cs 

There is a clear and unambiguous detection of the progenitor of 
SN2005cs in the ACS F814W image. As we have only one de- 
tection, the colour of the star makes the transformation from the 
F814W to standard Johnson I a little uncertain, and indeed all the 
transformations between the ACS filters and standard BVRI are 
subject to a similar uncertainty. However from the non-detections 
in the other bands, we can confidently say that it was a red star, later 
than mid K-type and the corrections between the ACS and Johnson 
VRI system are at most 0.2 mag, if the star was at the extreme 
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Figure 1. Pre- and post-explosion images of tlie site of SN 2005cs in M51. a) Post-explosion ACS/HRC F555W image (160 s) of SN 2005cs, the cross hairs 
indicate the centre of the PSF. The observation was split into short exposures to avoid saturation of the SN centroid. b) Pre-explosion ACSAVFC F555W 
image, in which the progenitor is not detected to mp555w = 25 at 5(T. c) Pre-explosion ACSAVFC F814W image, the progenitor is detected in this image, 
indicated by the cross hairs, with mFgi4w = 23.26. d) Pre-explosion NIRI J-band image, e) Pre-explosion NIRI //-band image, f) Pre-explosion NIRI 
A'-band image. The images are aligned to the reference ACSAVFC F555W image (panel b), such that North is up and East is left. 



M5 cool end ^Maund & Smartt'2005'). Hence in the rest of this dis- 
cussion we interpret the ACS filter magnitudes as standard BVRI, 
and the relatively small uncertainties are included in the error anal- 
ysis. The non-detection of a source in the BVR-band images is a 
very robust argument that the star was red. The spectral energy dis- 
tribution of the progenitor star in BVRIJHK is shown in Fig. 2. 
This is consistent with a red supergiant which is between K5 and 
M5. The star cannot have had bluer colours than a mid K-type, or 
it would have been detected in the V R-hwAs,. 

Although there is a range in possible spectral types that would 
fit the stellar SED, we can still set restrictive luminosity and mass 
li mits on the red prog enitor star. Using the standard reddening law 
of lCardelli eTall h989li) . it has a (V - /)o > 1.52, which c orre- 
spond s to a supergiant redder than approximately K3 l Elias et alj 
Il985h . The effective temperatures of red supergiants, ( V — I)o 
colours and bolome tric corrections from lElias et aljn985h and also 
iDrilling & Landoltl (2000 ) can then be used to determine a location 
box in the theoretical HR-diagram. This is shown in Fig. 3, with the 
stellar evolu tion tracks from the Geneva group plotted at Z — Zq 
iMevnetetal 111 994). The range in log L at a constant Tefj comes 
from the quadrature combined errors in distance, reddening, mea- 
sured 7-band magnitude (including colour correction), and bolo- 
metric correction (assuming plus or minus one spectral subtype). 
We can place quite a restrictive box on the location of the red su- 
pergiant progenitor in the HRD, even with a detection in only one 
band. The sta r was l ikely to have been a red supergiant of initial 
mass 9^2^10. |Liet_al. (2005c) report an initial mass of 7 — 9Mq, 
and indeed the position of the progenitor is closer to 7Mq in their 
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Figure 2. The spectral energy distribution of the progenitor is constrained 
by the /-band detection, and the upper limits in V RJ H K .The star can- 
not have been bluer than a mid K-type, or it would have been detected in 
the V and R b ands. The colours of the K5Ia-M4Iab Galactic supergiants 
were taken from' Elias et aljri985t) . and scaled to /=23.26 with appropriate 
reddening as discussed in the text. 



figure. The difference of 1 magnitude in the /-band is responsi- 
ble for the higher mass that we report, although within the errors 
the two studies are in agreement. This is a very similar mass to 
that derived for the progenitor of SN 2 003gd, which was also a 
red supergiant of initial mass SI^Mq jSmartt et alj Eooi. Both 
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Figure 3. The position of the progenitor of SN2005cs is restricted to He in 
the grey shaded region. The evolutionary tracks are the non-rotating models 
from the Geneva group (Meynet et al. 1994), where the initial mass, in units 
of Mq , is indicated for each track. The cool limit of the allowed region is 
the lowest observed T^g of late M-type supergiants. 



these supemovae were of type II-P \Ia et alJboOSd : iHendrv et all 
I2OO5I ) and future comparison of the photometric and spectral evo- 
lution of SN 2005cs will determine how similar th ey actually are. 
As pointed out by Smartt et al. ( 2004) and lLi et alj l2005cl. there is 
growing evidence that SNe II-P tend to arise in red supergiants with 
masses lower than ~15M0. There are restrictive mass limits, or de- 
tections, now set for seven SNe II-P (SNe I999br, I999em, I999gi, 
200 Idu, 2003gd, 2004et, 2005c s) and none of thein are higher than 
I5M(7) jM aund & Smartt 2005: ISmartt et all2003tlLi et all2005'3: 
IVan Dvk et al. 2003a b). One possibility often suggested to explain 
this is that the progenitors may be dust enshrouded red supergiants 
and hence appear to have lower mass and lower luminosity. Al- 
though the SNe themselves are not reddened, a plausible scenario 
is the destruction of this dust by the explosion jGraham & Meikld 
Il986h . However the deep i^-band image argues against this sce- 
nario. Even if the visual extinction were as high as Av ~ 5, the 
J-f— band images would be sensitive t o stars with Mk < —9.5, 
and with the bolometric corrections of lElias et alj Jl985h this im- 
plies log L/Lq < 4.6. Even when conservative errors on the dis- 
tance, limiting magnitudes and bolometric correction are included 
(resulting in ±0.2 dex). Fig. 3 shows that this rules out red super- 
giants with masses greater than I2M0. 

In summary, this letter presents the unambiguous detection of 
a red supergiant as the progenitor of SN 2005cs. We show that the 
likely initial mass for the star was in the range Mzams = 9^ . 
This adds to the emerging argument that all type II-P supemovae 
lie in the low to moderate mass range of red supergiants. The up- 
per limit from the deep Jf-band image is evidence that these are 
not dust enshrouded higher mass objects that have their cocoons 
destroyed in the explosion. 
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